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It has been suspected for nearly 50 years that clusters of galax- 
ies contain a population of intergalactic stars, ripped from galaxies 
during cluster formation or when the galaxies' orbits take them 
through the cluster center .§ | | Support for the existence of such 
a population of free-floating stars comes from measurements of 
the diffuse light in clusters Blilii, and from recent detections 
of planetary nebulae with positions and/or velocities far removed 
from any observed cluster galaxy. @0 But estimates for the mass 
of the diffuse population and its distribution relative to the galax- 



1 



ies are still highly uncertain. Here we report the direct detection 
of intergalactic stars in deep images of a blank field in the Virgo 
Cluster. The data suggest that approximately 10% of the stellar 
mass of the cluster is in intergalactic stars. We observe a rela- 
tively homogeneous distribution of stars, with evidence of a slight 
gradient toward M87. 

The process of cluster formation undoubtedly involved interactions and 
mergers of galaxies and tidal ablation of stars due to the mean gravita- 
tional shear of the cluster. E Further stripping is likely to have occured due 
to "galaxy harassment" , the cummulative effect of many close, high speed 
encounters with massive cluster galaxies. B Theoretical predictions for the 
fraction of mass in the intergalactic population range from 10% to 70%. This 
large uncertainty stems at a fundamental level from our poor knowledge of 
when galaxies formed, when and how clusters collapsed, and how the dark 
matter and baryonic matter were distributed within proto-galaxies. E3> a 
Detection and mapping of intergalactic stars could be an important tool for 
investigating these issues. 

Searches for the diffuse light from intergalactic stars are extremely diffi- 
cult due to the stringent requirements on background subtraction. Reported 
detections and upper lUwjijs range from 10% to 45% of the total cluster light 
in various clusters, a Q B a Even when a signal is detected, it is unclear how 
much is due to a true diffuse component, and how much to the extension of 
the galaxy luminosity function below the limits of individual source detec- 
tion. Intergalactic planetary nebula candidates have been identified in the 
Virgo cluster (3) and the Fornax cluster (10), but these numbers are not yet 
sufficient to draw strong conclusions on the distribution of intergalactic stars. 

In May 1996 we used the Hubble Space Telescope WFPC2 camera with 
the F814W (approximately I-band) filter to observe a blank field in the Virgo 
Cluster (Fig. 1). At faint limits compact galaxies are essentially indistin- 
guishable from stars, so the Hubble Deep Field, which was taken in an area 
of sky with no foreground cluster, was used as a control. For the com- 
parisons, the HDF images were divided into two groups, each consisting of 
exposures taken at three pointing positions with total exposure times iden- 
tical (to within 1%) to the Virgo Cluster field. Because the sky noise was 
considerably lower in the HDF images, additional noise was added to match 
the statistics in the Virgo Cluster field. These images were re-reduced and 
resampled to match the Virgo Cluster images. 

Figure 2 shows the comparison of source counts for the different images. 
Fainter than / = 26.5, there is a clear excess in source counts in the Virgo 
Cluster image. The excess amounts to approximately 630 extra sources to 
I = 27.9. We expect fewer than 20 Galactic foreground stars to this depth. 
The major uncertainty is the field-to-field variance in galaxy counts, because 
galaxies at the limiting depth of our image are only marginally resolved. For 
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a square field of angular size 6 on a side, the expected variance in the number 
counts is AN 2 = N + 2.2AN 2 A6 1 ~^ where N is the mean number per field, 
and the angular correlation function follows the powerlaw w{9) = AO 1-1 . EJ 
Using N from the HDF, and adopting A = 1.5 x 10~ 4 and 7 = 1.8, M the 
expected standard deviation in galaxy counts is AN = 80, well below the 
630 excess sources we detect. As a check that the HDF is not highly unusual, 
source detection and aperture photometry were carried out on an additional 
control field, centered on the radio galaxy 3C210. Source counts in this field 
match the HDF to within 20% down to the completeness limit (J = 27.5) of 
the 3C210 field. 

The total flux from the 630 excess sources detected in the image corre- 
sponds to a total magnitude / = 20.53, or a surface brightness (dividing by 
the area of the field) of \xj = 31.16 magnitudes per square arcsecond. The 
source density is relatively uniform across the field, with a marginally sig- 
nificant gradient towards M87. To estimate the surface brightness and total 
mass of stars below our detection limit, consider a population with a metal- 
licity, expressed as a decimal logarithm relative to thesolar iron-to-hydrogen 
ratio, [Fe/H] = -0.7, an age of 13 Gyr, a Salpeter |la initial mass function 
(N(M) oc M -2,35 ), arid a distance 18.2 Mpc. According to the theoretical 
luminosity functionO stars brighter than Mj = —3.4 (corresponding to a 
limiting magnitude I = 27.9 for the assumed distance) contribute 16% of 
the total flux from the stellar population. Our detection of 630 sources thus 
implies an underlying surface brightness of /i/ = 29.14. The theoretical color 
of the population is B — I — 1.97, so the _B-band surface brightness would 
be hb — 31.11; the surface mass density would be 0.14M Q pc if the initial 
mass function continues to O.1M . 

A purely empirical estimate of the surface brightness can be made by com- 
parison to previous HST observations of the galaxy NGC3379.H In those 
observations, approximately 39000 stars were detected in the three wide- 
field detectors, of which we estimate 5000 (±50%) were brighter than / = 
26.9 (roughly equivalent to our detection limit in Virgo; see below). The 
NGC 3379 surface brightness at the radius of these observations is \ib = 27.3. 
The difference in distance moduli between NGC3379 and the Virgo Cluster 
is 1.0 mag.EJ Thus our detection of 630 stars in Virgo implies a surface bright- 
ness — 2.51og(630/5000) + 1.0 magnitudes fainter, or \i B = 30.6, in reasonable 
agreement with the theoretical determination, given the uncertainties. 

The intergalactic stars are likely to have originated primarily from the 
elliptical and SO galaxies in the cluster. The early-type galaxies are more nu- 
merous in the cluster, have older stellar populations (and hence higher stellar 
mass-to-light ratios), and are likely to have inhabited the central megaparsec 
of the cluster for much longer than the spirals and irregulars.Ell To compare 
the intergalactic light to the light emerging from the early-type galaxies, we 
divide the Virgo cluster catalog Ea into successively larger annuli of 70 galax- 
ies centered on M87 and fit mean surface brightness vs. radius, where the 
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mean surface brightness is the total flux from galaxies divided by the area of 
each annulus. This fit gives hb = 28.7 at r = 44.5'. (This becomes /ib = 27.9 
if spirals and irregulars are included.) 

The implication is that about 10% of the stellar mass in the cluster is in 
the intergalactic component. 

While the distinction between intracluster stars and the M87 halo be- 
comes purely semantical at a radius of 235 kpc, it is of some interest to 
compare the inferred surface brightness to the extrapolated surface bright- 
ness of M87. An r 1//4 -law fit to the Caon et al. H photometry between 4.2 
and 9.2 arcminutes, predicts a surface brightness at r = 44.5 arcminutes of 
fiB — 30.3 along the major axis, and 33.4 along the minor axis. Our field is 
situated 25 degrees from the minor axis. Thus our detected number of stars 
is higher than expected from a pure extrapolation of the visible portion of 
M87. NGC4552, located 28.8 arcmin from our field, could also contribute 
to the source counts. Even if we assume no tidal truncation, a fit to the 
minor-axis photometry from 4.3 to 6.5 arcmin predicts = 32.5 at the po- 
sition of our field, which would correspond to ~ 25% of the detected source 
density. Overlapping halos of galaxies evidently account for only a portion of 
the stars we detect, unless the outer profiles deviate substantially from the 
inner r 1 / 4 laws. 

The estimated stellar mass fraction is sensitive to the age and metallicity 
of the population. For a distance modulus (M — m) = 31.3, metallicity 
[Fe/H] < —0.4 and age > lGyr, the underlying stellar mass fraction for 
our sample can plausibly range from 4 to 12%. At still younger ages the 
mass fraction goes down, while at still higher metallicities it goes up. Recent 
measurements of the Virgo Cluster distance modulus range from (m — M)q = 
31.0 to 31.7. tB The mass fraction is still ~ 5 — 10% unless the distance 
modulus is at the high end of this range and the metallicity is greater than 
[Fe/H] = -0.7. 

The present observations thus favor a rather smallmass fraction in the 
diffuse population. In contrast, Theuns and Warren Ej infer from the de- 
tection of 10 planetary nebula candidates, that 40% of the stellar mass is 
in the diffuse population of the Fornax Cluster. eJ Our results are consis- 
tent only if the intergalactic population is metal-rich, or the Virgo Cluster 
is more distant than ~ 20 Mpc. However, the planetary-nebula estimate 
suffers from small number statistics, possible contamination of the sample 
by background emission-line galaxies, and variations in the planetary-nebula 
specific frequency with stellar population. 

For an assumed distance of 18.2 Mpc, the total gravitating mass estimated 
from X-ray observations in the central 240 kpc of the Virgo Cluster is 1.2 — 
3.5xlO 13 M , E3 and the corresponding mass to light ratio M/ Lb ~ 260. Most 
of the light in the central region of the Virgo Cluster resides in the central 10' 
of M87, while most of the mass resides at larger radii. For the region outside 
this central 10', our observations imply a mean M/Lb ~ 700. The relatively 
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smooth distribution of mass inferred from the X-ray observations!! suggests 
that most of the intergalactic material was stripped via tidal interactions 
with the cluster potential, although we cannot rule out the possibility that 
some of the stars formed in situ, or that some were stripped off by impulsive 
interactions between galaxies ("harassment"). 

The luminosity of the tip of the red giant branch has been promoted as 
a distance indicator of comparable precision to Cepheid variables. @ While 
the number of stars in our sample is not large enough for a precise distance 
estimate, the shape of the observed luminosity function is consistent with 
that seen in the halo of NGC 5128 @, and in the halo of NGC 3379 @, shifted 
to a distance modulus in the range 31.2 < (m — M) < 31.6. 

Future observations should reveal whether the profile of the intergalac- 
tic stars is similar to that of the gravitating mass, whether it is centered 
on M87, and whether it is smoothly distributed. The new instruments on 
HST make it possible to measure the metallicity distribution of the diffuse 
stellar population and compare it to the metallicity distribution in the outer 
parts of elliptical galaxies and in the X-ray gas. Star counts also open the 
possibility of measuring the tidal radii of galaxies in the Virgo Cluster, and 
searching for the wakes of galaxies moving through the intracluster material. 
The combination of HST imaging and ground-based detection of planetary 
nebulae provides a powerful new probe of the cluster environment and of this 
nearly invisible population of stars. 
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Methods 

The zero-point of the instrumental magnitude scale was established by mea- 
suring the fluxes of the few bright isolated stais in the field in 0.5" radius 
apertures following the standard prescription, t^l Additional corrections were 
made as follows: (a) +0.05 mag to correct for the so-called "long exposure" 
effect (Hill et al. in press); (b) —0.04 mag, which is the appropriate color term 
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to convert to standard Cousin's I magnitudes for red stars with V— J~ 1.5;@ 
and (c) —0.04 mag to allow for the expected foreground extinction]^] We es- 
timate that the combined uncertainty from the zero-point calibration and 
the various correction terms amounts to 0.06 mag. Clearly extended sources 
were removed from the catalog by restricting our analysis to sources with a 
DAOPHOT sharpness parameter —0.6 < s < 0.4. A similar procedure was 
also followed for the HDF control field. Finally, we created 12 test data-sets 
in which 265 simulated stars were added to the real data frames. The input 
magnitudes of the simulated stars covered the range from 23.5 to nearly 30. 
We then processed these frames in the same way as the original data and 
produced a matrix relating the recovered magnitudes and detection efficiency 
to the input magnitudes of the simulated stars. From these simulations we 
estimate that the catalog is ~80% complete at / = 27.9, and use this as the 
limiting magnitude in our analysis. 
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Figure 1: 



The 'X' marks the position of the HST images in this view of the core 
of the Virgo Cluster. The solid circles mark positions of cluster galaxies, 
and are scaled by the square root of the galaxy luminosity. The central 
bright galaxy is M87, and the dashed circle marks a region 44.5 arcminutes 
from M87. The field is located on this circle ~ 235 kpc east of M87 (J2000 
coordinates 12 h 33 m 52 s , 12°21'40") and well separated from any of the other 
bright galaxies in the cluster (Figure [I]). The WFPC2 field covers 5.33 square 
arcminutes. The nearest galaxy to the WFPC field is a dE galaxy with 
a magnitude B=18.0, 6.6' away. The nearest bright galaxy is NGC4552 
(B = 10.8), at a distance of 23' (= 122 kpc if the distance to the Virgo 
Cluster is 18.2 Mpc). The total exposure time for the observation was 33500 
seconds. The exposures were taken at three positions separated by ±0.33" 
to improve resolution and aid in removing hot pixels and other cosmetic 
defects. Data were reduced and images combined following the procedures 
for the Hubble Deep Field (HDF). El The final reduced images were resampled 
via "drizzling" (Fruchter and Hook 1997, in preparation) onto a pixel grid of 
0.08". 
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Figure 2: 



Comparison of source counts in the Virgo Cluster field and the HDF. The 
histogram is the Virgo Cluster field. The points with error bars are the HDF. 
Magnitudes are based on PSF fitting. The excess in the Virgo Cluster field 
begins at / ~ 26.8. Source detection was performed using DAOFIND, and 
profile-fitting photometry was performed with the IRAF implementation of 
the DAOPHOT ALLSTAR code. See Methods for further details. 
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